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ABSTRACT 

We discuss the problem of using stellar kinematics of early-type galaxies to constrain 
the galaxies' orbital anisotropies and radial mass profiles. We demonstrate that com- 
pressing a galaxy's light distribution along the line of sight produces approximately 
the same signature in the line-of-sight velocity profiles as radial anisotropy. In partic- 
ular, fitting spherically symmetric dynamical models to apparently round, isotropic 
face-on flattened galaxies leads to a spurious bias towards radial orbits in the mod- 
els, especially if the galaxy has a weak face-on stellar disk. Such face-on stellar disks 
could plausibly be the cause of the radial anisotropy found in spherical models of 
intermediate luminosity clhpticals such as NGC 2434, NGC 3379 and NGC 6703. 

In the light of this result, we use simple dynamical models to constrain the outer 
mass profiles of a sample of 18 round, early-type galaxies. The galaxies follow a Tully- 
Fisher relation parallel to that for spiral galaxies, but fainter by at least 0.8 mag 
(/-band) for a given mass. The most luminous galaxies show clear evidence for the 
presence of a massive dark halo, but the case for dark haloes in fainter galaxies is 
more ambiguous. We discuss the observations that would be required to resolve this 
ambiguity. 

Key words: celestial mechanics - stellar dynamics - galaxies: kinematics and dynam- 
ics " galaxies: halos 



INTRODUCTION 



While it has long been known that spiral galaxies are embed- 
ded in halos of some kind of dark component, the evidence 
that ellipticals have dark halos is much more patchy. Anal- 
yses of the temperature profiles (e.g., Brighenti & Math- 
■ ews 1997), the mean temperature (Loewenstein & White 
1999) and the shapes (Buote & Canizares 1997) of the hot, 
X-ray emitting gas in bright ellipticals all demonstrate that 
these large galaxies have dark halos extending to many effec- 
tive radii. Similarly, results from gravitational lensing stud- 
ies (e.g., Kochanek 1996) show that the mass-to-light ratios 
of at least some large ellipticals must rise beyond an effective 
radius or so. It was unclear whether faint, compact ellipticals 
also had dark halos until the recent work of Rix et al. (1997) 
and Gerhard et al. (1998). These authors applied sophisti- 
cated spherical dynamical models to deep measurements of 
the line-of-sight velocity profiles (VPs) of the round galaxies 
NGC 2434 and NGC 6703. They found that the mass-to- 
light ratio in each rises outwards beyond about an effective 
radius. But despite this recent progress, our understanding 
of the properties of dark haloes of ellipticals remains poor. 

The main difficulty in determining the mass distribu- 
tions (or, equivalently, the potentials) of most ellipticals is 
the absence of dynamical tracers other than stars, plane- 
tary nebulae (PNe) and globular clusters. Since these trac- 
ers are coUisionless, one cannot determine the galaxy's po- 
tential without simultaneously constraining the distribution 
of orbits of the tracer population (the orbit distribution 



function, or DF) within the galaxy. There are two general 
approaches to this problem. One is to consider a range of 
parametrized trial galaxy potentials and, for each trial po- 
tential, to search for the DF that best matches the observed 
VPs. Potentials for which no satisfactory fit can be found 
can be ruled out. This is the approach used by Merritt & 
Saha (1993), Merritt (1993), Rix et al. (1997) and Gerhard 
et al. (1998). An alternative approach is to consider a range 
of simple parametrized forms for the DF and to search for 
the mass distributions consistent with the assumed DFs and 
the available observations. This has been used by Merritt & 
Tremblay (1993), Gebhardt & Fischer (1995) and Merritt 
(1996), but only considering one form for the DF (isotropy) 
and without making full use of the information contained 
in the observed VPs. We follow this second approach in the 
present paper, but using a wider range of parametrized DFs 
and making fuller use of observed VPs. 

The paper is organized as follows. In common with most 
other work in this area, our models assume spherical sym- 
metry for convenience. So, we begin in §2 by identifying the 
biases that this assumption introduces. In §3 we consider 
the simple problem of determining, for a spherical galaxy, 
the range of mass-to-light ratio profiles consistent with an 
observed velocity dispersion profile and some assumed form 
for the anisotropy of the velocity dispersion tensor. Section 4 
then shows how this method can naturally be extended to 
use VP shape information, including information from dis- 
crete tracers of the VPs, such as PNe. Applications to real 
galaxies follow in §5, before we sum up in §6. 
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2 THE DEGENERACY BETWEEN MASS, 
ANISOTROPY AND FLATTENING 

Throughout this paper we will assume that our round galax- 
ies are close to spherical. Before plunging into the details of 
the modelling procedure, let us first pause to consider the 
consequences of this assumption. It is instructive to take a 
real galaxy as an example. We use NGC 2434, an EO studied 
in detail by Rix et al. (1997). 

The points on Figure 1(a) plot the observed stellar kine- 
matics of this galaxy. In addition to the usual velocity dis- 
persion a, there is a Gauss-Hermite parameter /14 that mea- 
sures the deviation of the observed VP from the underly- 
ing best-fit Gaussian (van der Marel & Franx 1993). VPs 
more boxy than Gaussian have /14 < 0, whereas those more 
triangular have hi > 0. The dashed curves plot the kine- 
matics of a spherical model galaxy with the same projected 
light distribution as NGC 2434, an isotropic velocity distri- 
bution, and an isothermal potential with Vc = 300 km . 
The model reproduces the observed velocity dispersion quite 
well, but not the shapes of the VPs - its VPs are not cen- 
trally peaked enough. This can be rectified by increasing the 
fraction of stars on radial orbits in the model, while keeping 
the light distribution and potential fixed. Indeed, Rix ct al. 
(1997) have used sophisticated modelling machinery based 
on Schwarzschild's (1979) method to construct a strongly ra- 
dially anisotropic spherical model that provides an excellent 
fit to the observed kinematics of this galaxy. 

It turns out, however, that there is no need to invoke ra- 
dial anisotropy to explain these kinematics. Let us start with 
the spherical isotropic model again, but, instead of increas- 
ing the fraction of stars on radial orbits, let us keep the ve- 
locity distribution isotropic while flattening the model along 
the line of sight. Flattening the model like this increases the 
fraction of stars orbiting in the plane of the sky (i.e., orbits 
with high Lz) and decreases the model's projected velocity 
dispersion. This was first noted by Bender, Saglia & Ger- 
hard (1994), who invoked almost face-on circular orbits in 
a similar way to explain the minor-axis VPs of NGC 4660. 
The extra orbits in the plane of the sky affect the VP shapes 
in the same way as a bias towards radial orbits in a spherical 
galaxy. 

To give a more concrete example, we have increased 
the mass of our spherical isotropic model of NGC 2434 from 
Vc = 300kms~^ to 400kms~^, and squeezed its light distri- 
bution by first making it oblate with axis ratio q = 0.6 and 
then adjusting its cos 45 Fourier coefficients until its edge- 
on projected surface brightness distribution has the form 
shown in Figure 1(b). The 04 coefficient on this plot mea- 
sures the lowest-order symmetric deviation of each isophote 
from the underlying best-fit ellipse. The model is quite disky 
(04 > 0), with a peak 04/0 of 3.7% when viewed edge-on, or 
1.9% when averaged over all viewing angles. 

We calculate the kinematics of the model using the 
method described in Magorrian & Binney (1994). The re- 
sults are plotted as the heavy curves on Figure 1(a). This 
fiattened, isotropic model provides just as good a fit to the 
observations as the radially anisotropic, spherical models of 
Rix et al. (1997). While the diskiness of our model affects the 
shape of its VPs when seen face-on, even when seen edge on 
the disk does not have a strong effect on the model's (vo/cr), 
the characteristic ratio of the mean-streaming velocity to the 
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Figure 1(a). Kinematics of the EO galaxy NGC 2434. The points 
on plot the observed velocity dispersion a and Gauss-Hermite 
VP shape parameter hi versus radius. The broken curves plot 
the corresponding quantities for spherical isotropic galaxy mod- 
els. The heavy solid curves show the kinematics of an flattened 
isotropic model of the galaxy, viewed face-on. Figure 1(b) shows 
the edge-on surface brightness distribution of this model. 

velocity dispersion (Binney 1978; Davies et al. 1983). This 
remains very close to the predictions for an oblate isotropic 
rotator. 

Of course, these results do not mean that NGC 2434 
must be isotropic with a weak face-on stellar disk, but they 
do show that making unwarranted assumptions about a 
galaxy's intrinsic shape will lead to biases in the orbit dis- 
tribution that one finds. Throughout the rest of this paper 
we will model our round galaxies under the assumption of 
spherical symmetry. In Appendix A we quantify the degen- 
eracy between mass, anisotropy and shape, showing how the 
results from our spherical models may be corrected for any 
supposed degree of (constant) flattening. 

3 CONSTRAINING THE MASS PROFILE 
USING VELOCITY DISPERSION 
MEASUREMENTS 

We are given measurements of a galaxy's projected velocity 
dispersion ai with corresponding errors (assumed Gaussian) 
AiTi at points J?,, with i = 1 . . . N . The luminosity density of 
the galaxy is j{r). We seek the range of smooth mass-to-light 
ratio profiles T(r) consistent with these observations and 
some assumed form for /3(r), the anisotropy of the galaxy's 
velocity dispersion tensor, where 

and (Jr(r) and <J<j,{r) are the velocity dispersions in the ra- 
dial and azimuthal directions at radius r. Before explaining 
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Figure 1(b). Edgc-on surface brightness distribution of the flat- 
tened isotropic model in Fig. 1(a). The panels show the model's 
isophotal shape parameters (surface brightness, ellipticity and 
lowest-order shape parameter 04/0), using the convention of Ben- 
der & Mollenhoff (1987). 

how we determine the range of T(r), we first describe the 
necessary simpler step in the method: how to calculate the 
predicted velocity dispersion profile &{R) of a trial T(r). 

3.1 Prom T(r) to a{R) 

Given T(r), the mass interior to radius r is given by 



M{r)=4n T{r')j{r')r'^ dr' . 
Jo 



The Jeans equation (e.g., Binney & Trcrnaine 1987), 



d{jV: 



dr r 



.GM 



(2) 



(3) 



can be integrated to give the intrinsic second-order moments 
in terms of j{r), M[r) and /3(r): 



J r 



jvf{r)= I drv(r) — ^ — exp 



2/3(r"; 



dr ' 



and, by definition. 



vl = vl = {l-P)vl 



(4) 



(5) 



Projecting along lines of sight, the zeroth and second-order 
moments at projected radius R are 



I{R) = 2 



j{r)r dr 
rdr 



n2 

J?2 V ' r2 



(6) 



To include the efi^ects of seeing, we assume that the PSF 
of the kinematic observations is a Gaussian with dispersion 
fT* = FWHM/2.305, 



PSF(Aa;) 



27ra2 



exp 



1 /AxV 
'2 \~^) 



(7) 



The convolution of any projected distribution /Lio(i?) (e.g.. 
I{R) or /wf(7?)) with this PSF is 



Jo 



R exp ( -—^ 



lo 



^ ] ^o{R')dR', 



(8) 



where Io{x) is the zeroth-order modified Bessel function. Di- 
viding the seeing-convolved second-order moment by the ze- 
roth order one and taking the square root yields the model's 
predicted cr(R), and, in particular, its predictions ai at the 
observed points Ri. 

In the numerical implementation of the calculations 
above we store the logarithms of j{r), M{r), etc. on grids 
spaced logarithmically in radius from the radius of the in- 
nermost observed isophote to a few times the radius of the 
outermost one, typically with about 30 grid points. Values 
of j{r), M{r), etc. at intermediate points are obtained by 
linear interpolation on these log-log grids and exponentiat- 
ing. Since the positions of the radial grid points, as well as 
j{r) and /3(r), are independent of T(r), many of the coeffi- 
cients involved in the numerical calculation of the integrals 
(2), (4), (6) and (8) need only be evaluated once, leading to 
a considerable decrease in the time needed to calculate the 
ai from a trial T(r). 

3.2 Determining the range of acceptable T(r) 

By Bayes' theorem, the probability of any model T(r) given 
the available velocity dispersion data D and our assumed (3 
satisfies 



v = pir\D,i3)cxp{D\r,i3)pir). 



(9) 



Since we assume that the observational errors are Gaussian, 
the likelihood p{D\T,f3) is given by 



p(i;|T,/3)=exp (-ix 



where 



N 
i=l 



(10) 



(11) 



The prior p(T) encodes our prejudices about what consti- 
tutes a reasonable T(r). Our only expectation is that it be 
smooth, and we take 

Inp(T) = (^^) dlnr • J dlnr, (12) 

which uses the mean-square change in the slope of In T ver- 
sus In r as a measure of smoothness. For a reasonable T(r) 
the change in the mean-square slope would be less than 
about 100. We set A — 1/100 so that an increase in the 
mean-square slope by 100 has the same effect on p(T|D,/3) 
as an increase in of one. 
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To obtain the range of T(r) that maximize 'P, we first 
fit an initial guess of the form 

T« = To(^)"(l+^)\ (13) 

using the downhill simplex method (Press et al. 1992) to 
find the parameters (To,ro,a, b) that minimize in equa- 
tion (11). Then the Metropolis algorithm (Metropolis et al. 
1953) is used to explore the range of T(r) consistent with 

the data: 

(1) Calculate V\T] for the initial guess. 

(2) Make a change <5T to T. Calculate x^[T + 5T], p(T + 
5T) and thus7'[T-|-5T]. 

(3) if AP = P[T + ST] - ^[T] > accept the change 5T- 
otherwise accept it with probability exp(A'P). 

(4) Go back to step 2. 

Each change 5T in step (2) is made by choosing a sin- 
gle radial grid point i at random and adding xAi to InT,, 
where x is a random number uniformly distributed in (—1, 1) 
and the constant A, sets the size of the maximum change in 
InT at the point i. For maximum efficiency the A, should 
be chosen so that the change in each step of the Metropo- 
lis algorithm has roughly equal chances of being accepted 
or rejected. We ensure this is the case by initially choosing 
all A, = 1 and following the Metropolis algorithm from the 
initial fit (13) for enough iterations to ensure that each ra^ 
dial grid point is sampled many times. During this time, we 
increase A^ by a small factor (e.g., 1.5) if a change in the 
j*h grid point is accepted, and decrease it by the same fac- 
tor if the change is rejected. We then restart the Metropolis 
algorithm from the initial fit (13) keeping the A, fixed. 

As the Metropolis algorithm progresses, the probabil- 
ity of any given configuration T(r) occurring tends towards 
the probability p(T|_D,/3) = V ■ This makes it easy to ob- 
tain confidence intervals on T(r) and other model quanti- 
ties. In addition, if we have two models, one with anisotropy 
/3i(r), the other with /32(r), and no reason for favouring one 
over the other (p(/3i) = p(/32)), then a simple application of 
Bayes' theorem yields their relative probability, 

p(/3i|g) ^ P(g|/3i) ^ /p(g|T,/3i)p(T)dT 
v{p2\D) p(Z)|/32) /p(Z)|T,/32)p(T)dT' ^ ^ 

which from equation (9) is just the ratio of the average values 
of V in the Metropolis algorithm. Of course, equation (14) 
is not very useful in the present case, since for most choices 
of (3{r) there is usually some smooth T(r) that projects to 
a good fit to the observed CTj. But given further information 
on the shapes of the VPs (Section 4), and if we somehow 
knew that the galaxy was close to spherical, we could use 
equation (14) to distinguish between models with different 
Pir). 

3.3 An application to M87 

In principle, the mass profile can also be obtained by rear- 
ranging the Jeans equation (3) to get 

in which Vr follows from an inversion of equation (6). This 
approach was used by Sargent et al. (1978) to obtain the 



mass distribution of M87 assuming that it is isotropic (/3 = 
0). The method is very direct, but requires the fitting of a 
smooth, arbitrary a{R) profile to the observed Ui, leading 
to complex biases in the results. Here we use Sargent et al.'s 
data and results to provide a simple test of our method. 

We obtain our model luminosity density j{r) from 
their photometry using the method described in Magorrian 
(1999). It is unclear exactly what seeing was in effect at 
the time of Sargent et al.'s observations, but wo assume 
2 arcsec FWHM. We then apply the algorithm above to 
their kinematical data for 10'' iterations, the results of which 
are plotted on Fig. 2. The panels show the 95% confidence 
limits on the local mass-to-light ratio T(r), the "average" 
mass-to-liglit ratio M{r)/ L{r) enclosed within radius r, and 
the corresponding limits on the model's projected dispersion 
profile a(r), together with the observed a,. The values of 
M(r)/L(r) found by Sargent et al. are also plotted on Fig. 2 
for comparison with ours. There is quite good agreement 
between our results and theirs, except within the innermost 
few arcsec of the galaxy centre. The reason for this is simple 
(e.g., Binney & Tremaine 1987): at large radii, the j{r) pro- 
file is steep (~ r~^) while the profile is shallow (~ r°), so 
that the value of the bracketed expression in (15) used by 
Sargent et al. is dominated by the relatively well-determined 
slope of j{r). At smaller radii, j{r) becomes much shallower, 
Vr steepens, and so Af(r) is affected more by the poorly de- 
termined slope of the latter. 

While the agreement between the results from the two 
methods is gratifying, it is important to note the fundamen- 
tal difference between them. In methods like Sargent et al.'s 
(e.g., Gebhardt & Fischer 1995), one is required to make as- 
sumptions about unobserved data, whereas in methods like 
ours (e.g., Merritt & Tremblay 1993; Merritt 1996) one uses 
only real, observed data, and confines any assumptions to 
the model itself. In the former case, there is no straightfor- 
ward way to obtain confidence bounds on the model's mass 
density profiles, or indeed to impose the reasonable con- 
straint that the mass density must be non-negative. On the 
other hand, our model's disregard for imaginary data is why 
the exact amount of mass enclosed within the innermost few 
arcsec is quite uncertain - due to the effects of seeing, the 
data are consistent with quite a gentle increase in mass-to- 
light ratio near the centre. Similarly, our model correctly 
reproduces the fact that there are no interesting constraints 
on T(r) beyond the radius of the outermost observed dis- 
persion. Of course, the one big advantage of methods like 
Sargent et al.'s over ours is that they are relatively quick, 
and can be adequate if one is content with finding only a 
restricted range of the full set of plausible mass density pro- 
files. 



4 USING VELOCITY PROFILES TO 
CONSTRAIN BOTH THE MASS PROFILE AND 
THE ANISOTROPY 

The method above can be extended to make use of the full 
observed VPs if, instead of assuming a form for 0{r), we 
make a stronger assumption about the form of the galaxy's 
distribution function (DF). By Jeans' theorem, the DF of 
a spherical galaxy is a function only of the stars' binding 
energies £ and angular momenta L per unit mass. We make 
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Figure 2. Our model for the mass distribution of M87, using the 
data of Sargent ct al. (1978) and assuming the galaxy is isotropie 
(/3 = 0). The bottom panel shows the 95% confidence limits on 
the the local mass-to-light ratio T(r). The corresponding limits 
on the "global" mass-to-light ratio M(r)/L(r) and the model's 
projected velocity dispersions are plotted on the pajiels above. For 
comparison, the points show the M(r)/L{r) obtained by Sargent 
et al. (middle panel) and the observed dispersions (top panel). 



the following ansatz for the DF (Cuddeford 1991): 



f{£,L^) = L-^^<'fo{Q) withQ = f--^. (16) 

We assume that fo{Q) = for Q < 0. The luminosity den- 
sity of these models 



r(f-/3o) 



Jo ■ 



MQK^ - Q)'^'-^" dQ, 



(17) 

where ^(r) is related to the gravitational potential <I>(r) 
through 4>{r) = — $(r). This demonstrates that the func- 
tion fo{Q) is (implicitly) determined once the luminosity 
density j{r) and mass distribution are specified. 

More generally, the non-trivial velocity moments of 
these models are given by 



[v^r''vfvlmr)=jv'r''vfvf-{r) = Cklmir)Vk+l+m{r), (18) 



v„ir)= Vo(g)(V'-Q) 

Jo 



ri-/3o-H/2 



dQ, (19) 



and the coefficients 



2X -(l+m+l-po) 



B 



hi 



l + m + l 



(20) 

This helps to clarify the moaning of the parameters /3 and rai 
from (18) and (20), the anisotropy 



af + ri 

which tends to /3o as r — » 0, and to 1 for r » r^- 



(21) 



4.1 From T(r) to VPs 

It is straightforward to extend the method of section 3 to 
calculate the higher-order moments (18) and reconstruct the 
VPs from them. Our method avoids the difficult problem of 
obtaining fo{Q) explicitly from (17). It closely follows the 
approach of Magorrian & Binney (1994). 
Integrating (19) by parts gives 



Vnir) = (n-/3o + l/2) / Vr 



-t(r')^d/. 



(22) 



which is the higher-order analogue of the Jeans' equa- 
tion (3). Having Vo from j{r), it allows the full set of V„ 

up to any arbitrary order to be calculated. The moments 
[v^''vfvj"']{r) then follow from (18) and (20). Projecting 
the n"^-order moments along lines of sight and convolving 
with seeing as in Section 2 yields the classical moments of 
the VPs C{R;Vp), namely 



[v;]{R) = Iv^{R) = I{R) f C{R;Vp)v;dvp. (23) 

J — oo 

Each VP can be reconstructed from its moments Vp using a 
Gram-Charlier series of type A (Kendall & Stuart 1943): 



£(^p) = M!f^V*/^eiK), 
(7e ^-^ 



(24) 



where Wg = Up/ug, a{x) = exp{—x'^/2)/^/2Tv is the standard 
Gaussian, the Hei are Hermite Polynomials of the first kind 
(e.g., Abramowitz & Stegun 1965), and the coefficients 



£{Vp)Hei{Wg)dVp 



(25) 



are just linear combinations of the Vp. Kendall & Stuart 
(1943) give conditions sufficient for convergence. For realistic 
VPs to be well-approximated using a small number of di, the 
expansion variable ag should be chosen so that ag ~ 

As in §3.1, many of the coefficients involved in the nu- 
merical ovalution of the V„, Vp and di depend only on /Jo, ra 
and j{r) and therefore need only be calculated once, signif- 
icantly speeding up the calculation of the di{R) from T(r). 
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4.2 Comparing model and observed VPs 

Since VPs are usually close to Gaussian, observed VPs are 
usually reported by giving the first few terms of their Gauss- 
Hermite expansion (van der Marel & Pranx 1993; see also 
Gerhard 1993): 

C{v^) = WHl y h,iH2iiw,), (26) 

where Wg = Vp/ag, and the Hi arc Hcrrriite polynomials of 
the second kind. There is a different set of Gauss-Hcrmite 
coefficients (/lo, ^2, /i4, • • •) for each choice of (7g,crg), but 
observers fix the latter to be the line-strength and disper- 
sion of the best-fitting Gaussian to the VP. Then ho = 1, 
h2 = and the lowest-order deviation of the VP from the 
best-fitting Gaussian is measured by h4 (van der Marel & 
Pranx 1993). 

While the most straightforward course of action would 
be to compare models and observations using this choice of 
(7g,crg), this has the minor disadvantage that it can only 
give model predictions hi at the same radii as the obser- 
vations. To remedy this, we instead choose 7g = 1, and 
CTg(i?) equal to the smoothly varying dispersion of the best- 
fit initial model (13). We use the relations (Magorrian & 
Binney 1994) 

dhi 1 , 

— = 

07g 7g 

1^ = ^ (v'(' + 2)(/ + l)fti+2 +hi- - l)hi-^ , 

(27) 

to make the small adjustments required to bring the Gauss- 
Hermite coefficients of the observed VPs to this choice of 
(7g, crg(Ji)). Conversely, when presenting our final results, 
we will use these relations to return the model VPs in the 
observational paraiiietrization (7g, o"g, /14, he, . . .) by adjust- 
ing (7g,a-g) so that (/io,ft2) = (1,0). 

The Gauss-Hermite coefficients of the model galaxy, 
hi{R), can be obtained from the moments Vp using the 
Gram-Charlier expansion (24) and some judicious manip- 
ulation of integrals of Hcrmite polynomials. Magorrian & 
Binney (1994) show that there is a simple linear relation 
between the two, namely 

00 

hk='^Ckidi, (28) 
1=0 

and describe how to obtain the matrix Ckii'jg, ag). Of course, 
in practice it is neither possible nor necessary to include all 
the terms in these infinite sums. Limited spectral resolution 
together with the effects of template mismatch mean that 
observations can only yield sensible information on the hi 
up to i = Nh — 4. For most of our models we find that 
the error in the predicted /14 introduced by truncating the 
series (28) at Z = 10 is of order 0.002, much less than typ- 
ical observational uncertainties. As a stronger test of our 
moment-based VP calculations, we have calculated the VPs 
of various models by solving eq. (17) numerically and inte- 
grating directly over the DP. The /i4 coefficients calculated 
this way differ from the moment-based results by less than 
0.002. 

To first order, errors in the hi are uncorrelated (van 
der Marel & Pranx 1993). Assuming that these errors are 



Gaussian, the likelihood of the trial T(r) is given by 

p(D|T,/3o,ra) = exp(-ixl) (29) 

with 

xS[T,ft,..i = t£(MMz|iM)'. (30) 

i=l j=0 ^ ' / 

In practice we truncate this series at ji — 2, since the un- 
known flattening of the galaxy affects its /1.4 profile (i.e., 
we calculate the full VPs, but fit only to a). Our method 
can easily be modified to deal with other parametrizations 
of VPs, including so-called "non-parametric" parametriza- 
tions. In most cases, however, a proper treatment of the 
correlations between the parameters is likely to lead to a 
significant increase in the complexity of the expression for 
the likelihood. 

4.3 Including discrete tracers of the VPs 

The rapid fall-off in the surface brightness of galaxies means 
that it is very difficult to use stellar absorption lines to mea- 
sure VP shapes reliably far from the galaxy centre. Fortu- 
nately, in any galaxy there will be many stars that have 
evolved into planetary nebulae (PNc). These objects emit a 
significant fraction of their light in the O [III] line at 5007A, 
making it easy to measure their line-of-sight velocities. If we 
make the reasonable assumption that PNe axe drawn from 
the same DF as the stars, then the probability that a PN at 
projected radius R will have line-of-sight velocity between Vp 
and Vp + dvp is just C{R; Vp) dvp. Therefore the likelihood of 
measuring A'pn PNe, the i*^ at radius Ri with line-of-sight 
velocity Vi and error Avi, 

iVpN 

p(Z)pn|T, 00, ra) oc C{Ri; Vi)Avi, (31) 

i = l 

assuming that the observational uncertainties Avi <C cr{Ri). 
Then, using the results from §4.2 above, the combined like- 
lihood of the PNe and VP data is 

p{D\r, /3o, ra) = p(L'vp|T, po, ra) p(L'pn|T, /3o, ra) 

« exp {^--^XhJ JJ jC{Ri,Vi)Avi. 

Since the shape of C{R,v) is affected by the galaxy's un- 
known flattening, in practice we simply ignore the deviations 
of the model's predicted C{R,v) from their best-fit Gaus- 
sians, consistent with our treatment of ft.4 measurements 
above. 

4.4 Summary 

In summary then, the procedure we use to obtain the range 

of T(r) consistent with observed VPs is as follows: 

(1) Fit an initial guess of the form (13) to the observed 
dispersions di. 

(2) Use (27) to adjust the observed Gauss-Hermite coeffi- 
cients hi to 7g = 1 and Og{R) equal to the dispersion 
profile of this smooth initial fit. 

(3) Iterate the Metropolis algorithm as in §3.2, but includ- 
ing higher-order (typically up to 10*''-order) moments 
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(§4.1) rather than just the second-order moments (§3.1). 
Use the Ukehhood (32) instead of (10). 
(4) When the Metropolis algorithm is finished, calculate 
the confidence limits on the model T(r), and the cor- 
responding hi. Convert the latter back to the observa- 
tional paranictrization (7g, a^,hi, . . .) using (27). 
If we somehow knew that the galaxy was close to spherical, 
we could constrain its anisotropy by constructing models 
with a variety of Ba and and using equation (14) to find 
the range of (/^o, Ta) preferred by the available data. In prac- 
tice, however, wc construct models for just a few choices of 
(/So, ''a), compare the predicted a and Ha profiles with the ob- 
servations and use the results of Appendix A to estimate the 
allowed range of flattening or anisotropy in the real galaxy. 

5 APPLICATION TO REAL GALAXIES 

In this section we apply the machinery to a number of real 
galaxies, listed in Table 1. The sample consists of galaxies no 
flatter than E2 on the sky, with published kinematics mea- 
sured out to of order the galaxy's effective radius or further. 
For most of the galaxies the available data consist simply of 
a velocity dispersion profile, but in a handful of cases mea- 
surements of the shapes of the VPs (hi parameters) are also 
available. Two of the galaxies have PNe radial- velocity mea- 
surements: 15 radial velocities for M32, 29 for NGC 3379. 

We obtain the luminosity density j(r) of each galaxy by 
deprojecting its spherically averaged surface brightness dis- 
tribution using the method described in Magorrian (1999) 
with rir = 30 radial grid points. Then we apply the ma- 
chinery of §4 using the observed kinematics to constrain the 
galaxy's mass profile, first under the assumption that the 
galaxy is isotropic (/3 = 0), then assuming that it is mildly 
radially anisotropic (/3 = 0.3). In each case we carry out 10 ' 
iterations of the Metropolis algorithm, including terms up 
to Z = 10 in equation (24). But, given the degeneracy be- 
tween anisotropy and flattening identified above, we ignore 
any measurements of VP shapes, and fit simply to the ve- 
locity dispersion profile (i.e., wo truncate the series (30) at 
Nh = 2), and merely predict the hi profile for a spherical 
galaxy with the assumed anisotropy. 

5.1 Results for individual galaxies 

The results of this modelling are plotted on Figure 3. One of 

the most obvious features of the results is that each model's 
outer mass profile is almost independent of the anisotropy 
we assume (Efstathiou, EUis & Carter 1980; Gerhard 1993). 
Some of the galaxies deserve more detailed comment. 

M32 Our models use the minor-axis velocity disper- 
sion profile obtained by van der Marel et al. (1994) and 
the 15 PNe radial velocities measured by Nolthenius & 
Ford (1986). These data are insufficient to place any in- 
teresting constraints on whether M32 has a halo or not: 
the Vc{r) and T(r) profiles of our models are both con- 
sistent with being fiat. The results plotted are for models 
where we have subtracted a heliocentric systemic velocity 
of -203 kms"^ (TuUy 1988) from the PNe velocities. As- 
suming systemic velocities of — 195kms~^ or — 185kms~^ 
(Nolthenius & Ford 1986) yields almost identical results. 

NGC 2434 As mentioned in §2, Rix et al. (1997) 
have used Scharzschild's method to construct spherical mod- 
els of this galaxy. Their models required a dark halo with 



Vc ^ 300kms~^ at around reS- Our results are in very good 
agreement. We re-iterate, however, that if this galaxy were 
flattened along the line of sight these spherical models would 
underestimate the true circular velocity. 

NGC 3379 Our models use the minor-axis dispersion 
profile measured by Statler & Smecker-Hane (1999), and the 
29 PNe radial velocities from Ciardullo, Jacoby & Dejonghe 

(1993) . We have subtracted a systemic velocity of 900 km s~^ 
from the PNe velocities for the model plotted in Fig. 3. The 
results are almost identical if we instead assume a systemic 
velocity of 916kms~^, equal to the mean recession veloc- 
ity of the PNe sample, or if we use the 881kms~^ quoted 
in TuUy (1988). Statler & Smecker-Hane (1999) discovered 
that this galaxy has a kink in its dispersion profile at a ra- 
dius of ~ 15 arcsec. Our models fit the data beyond the 
kink by having T(r) increase by a factor of three between 
15 arcsec and 60 arcsec, then becoming further out in order 
to fit the PNe radial velocities. An alternative explanation is 
that NGC 3379 may be a face-on SO (e.g., Capaccioli et al. 
1991), or a weakly tricixial system viewed face on (Statler & 
Smecker-Hane 1999). We note that the the evidence for posi- 
tive hi in this galaxy could be due either to radial anisotropy 
(Gebhardt et al. 2000) or to the galaxy being fiattencd sys- 
tem viewed face on (Section 2). 

NGC 4472 Brighenti & Mathews' (1997) modelling 
of the X-ray emission from this galaxy indicates a mass of 
about 3x 10^^ Mq within 7.8 kpc (scaling to our assumed dis- 
tance of 15.3 Mpc), which corresponds to Vc — 410kms~^. 
Similar results are found by Irwin & Sarazin (1996). This 
is in good agreement with our Vc = (419 ± ll)kms~'^ at 
Ireff = 7.8 kpc. 

NGC 4486 Nulsen & Bohringer (1995) find that the 
X-ray data for M87 are consistent with a mass per unit 
length of 3.6 x 10^° kpc"S or Vc = 390 km s"^ Similarly, 
Tsai (1993) finds a mass 2.8 x lO^^A/g enclosed within a ra- 
dius 7.8 kpc, also giving Vc = 390 km s^^. Using the velocity 
dispersion measurements of Sembach & Tonry (1996), our 
models find Vc = (530 ± 9)kms~'^, larger by 30%. Part of 
this discrepancy is due to the fact that Sembach & Tonry's 
measurements are systematically offset by 10 to 15% (Sem- 
bach, private communication). Comparison with the mea- 
surements of Sargent et al. (1978), Bender, Saglia & Gerhard 

(1994) and van der Marel (1994) confirms that Sembach & 
Tonry's measurements are systematically overestimated by 
about 15%. 

Taking our models at face value, five of the galaxies in the 
sample (M87, NGC 2434, 4472, 5846 and 7796) show evi- 
dence for a dark halo (flat rotation curve and rising mass-to- 
light ratio at large radii). Five (NGC 1379, 4434, 4552, 7192, 
7507) show evidence against a dark halo (constant inass-to- 
light ratio and falling rotation curve). It should be borne 
in mind, however, that these results depend critically on 
our assumption of spherical symmetry. For example, in our 
model of NGC 4552 the rotation curve falls off at large radii, 
and mass follows light very closely. Despite being round on 
the sky, however, this galaxy is classified as an SO in some 
catalogues. If its light distribution did indeed become signif- 
icantly flattened at large radii, then it would have to have a 
rising mass-to-light ratio proflle, and possibly a flat rotation 
curve. 

The remaining eight galaxies (M32, NGC 1400, 3379, 
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Table 1. Sources of data used. 



Name 


Dist. 




log Lb 


B - V 


V - I 




Phot. 


Kin. 


^max 






(Mpc) 










(arcsec) 


Source 


Source 


(arcsec) 


/I — 1 \ 

( km s ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


M32 


0.8 


8.76 


8.49 


0.89 


1.14 


40 


K87 


vdM94 


22 


75± 15 


















NF86 (PNe) 


271 




NGC 1379 


17.9 


11.66 


10.03 


0.91 


1.19 


24 


F1H89 


GCBZL98 


61 


167 ±5 


NGC 1400 


21.5 


11.62 


10.21 


1.03 


1.37 


38 


L85 


BBBDDSSZZ94 


45 


268 ± 20 


NGC 2434 


20.0 


11.43 


10.22 


0.91 


1.18 


24 


CD94 


CZMDQ95 


62 


296 ± 17 


NGC 3379 


9.9 


10.43 


10.01 


0.99 


1.23 


35 


PDIDC90 


SS99 


78 


262 ± 14 


















CJD93 (PNe) 


209 




NGC 4434 


15.3 


12.83 


9.43 


- 


- 


19 


CCR90 


SP97 


17 


142 ± 18 


NGC 4161 


15.3 


13.61 


9.12 


0.89 


1.19 


6 


CCR90 


SF97 


17 


165 ± 13 


NGC 4167 


15.3 


14.81 


8.64 


- 


- 


10 


- 


BN90 


9 


98± 10 


NGC 4472 


15.3 


9.32 


10.83 


0.99 


1.26 


105 


PDIDC90 


FIF95* 


140 


419 ± 23 


NGC 4478 


15.3 


12.14 


9.71 


0.90 


1.19 


14 


CCR90 


DEFIS83 


27 


220 ± 13 


M87 


15.3 


9.52 


10.75 


0.98 


1.27 


105 


PDIDC90 


ST96 


170 


531 ± 17 


NGC 4494 


15.3 


10.69 


10.29 






45 


BM87 


JS89* 


37 


240 ± 25 


NGC 4552 


15.3 


10.84 


10.23 


0.97 


1.20 


30 


MM94 


SP97 


58 


330 ± 17 


NGC 5846 


26.0 


10.67 


10.75 






63 


MM94 


ST96 


145 


418 ± 20 


NGC 6407 


60.0 


12.90 


10.59 


1.12 


1.38 


33 


CD94 


CD94* 


43 


440 ± 32 


NGC 7192 


34.0 


12.04 


10.44 


0.97 


1.24 


28 


CD94 


CD94* 


43 


308 ± 23 


NGC 7507 


21.0 


11.15 


10.38 


1.00 


1.28 


31 


F1H89 


BBBDDSSZZ94 


66 


321 ± 16 


NGC 7796 


37.0 


12.28 


10.42 


1.00 


1.23 


37 


SWJTB91 


BBBDDSSZZ94 


37 


371 ± 30 



Distances (col. 2) arc taken (in decreasing order of preference) from Faber et al. (1997), Faber ct al. (1989) or the Hubble relation 
assuming Hq = 80kms~'^ Mpc~^. B magnitudes, B — V colours and effective radii (cols. 3, 5 and 7) are taken from Faber et al. (1989). 
V — I colours (col. 6) are from Prugniel et al. (1993) and Poulain & Nieto (1994). The photometry used (col. 8) is supplemented with 
HST photometry from Faber et al. (1997) when available. Col. 10 gives the radius of the outermost observed point of the kinematic 
source (col. 9). Our models use velocity dispersions along one axis (usually the minor if available), except for sources marked with 
an asterisk for which we use results from multiple slit positions. Col. (11) gives the 95% confidence intervals on the circular velocities 
of the isotropic models at one effective radius. 



4464, 4467, 4478, 4494 and 6407) are indeterminate, with 
both a fiat mass-to-light ratio and a flat rotation curve at 
the outermost available radius. The most sensible way of 
distinguishing between constant-T and constant-Vc models 
of these galaxies would be by using PNo to estimate the 
galaxy's mean luminosity-weighted velocity dispersion out- 
side the radius of the outermost stellar velocity dispersion. In 
Appendix B we show that distinguishing between constamt- 
T and constant-Vc models with 95% confidence would re- 
quire radial velocities of ~ 100 such PNe. 

5.2 The Tully— Fisher relation for ellipticals 

Figure 4 plots the circular velocity derived by our models 
at Ireff versus galaxy luminosity. They lie on a Tully-Fisher 
relation parallel to that for spiral galaxies, as was first noted 
by Pranx (1993). Our elliptical Tully-Fisher relation is off- 
set from Pierce & TuUy's (1992) spiral relation in the sense 
that, for given Wc, ellipticals are fainter than spirals by 1.5 
mag {B band) or 1.0 mag (/ band). The offset from Mathew- 
son, Ford & Buchhorn's (1992) sample is 0.83 mag in I. Of 
course, these numbers are underestimates (overestimates) of 
the true offset if our ellipticals are oblate (prolate) galaxies 
viewed face on. 

The scatter in our Tully-Fisher relation for ellipticals 
is 0.55 mag, significantly larger than the values ~ 0.2 mag 
reported for samples of nearby spirals. But the scatter in 
the real v^-L correlation for ellipticals could well be smaller 
than 0.55 mag: our galaxy sample is heterogenous, preclud- 
ing accurate distance (and therefore luminosity) estimates. 



and we are unable to correct our circular velocities for the 
unknown flattening of our galaxies, a problem that is usually 
unavoidable when using elliptical galaxies. 

Our results differ from those of Neistein et al. (1999), 
who looked for a Tully Fisher relation in a sample of nearby 
SO galaxies, using the method of surface brightness fluctua- 
tions to obtain careful distance estimates, and obtaining cir- 
cular velocities by deprojecting the galaxies' observed stellar 
kinematics. They found a relation with an intrinsic scatter 
of 0.7 mag (/ band), with SOs of given Vc only 0.5 mag 
fainter than corresponding type-Sc spirals. This is less than 
expected if SOs are former spirals with truncated star for- 
mation histories. Our results can only be brought into broad 
agreement with theirs if we assume that our round ellipticals 
are mostly face-on prolate systems. Their method, however, 
also requires assumptions about the intrinsic shapes of their 
galaxies, although it is unclear whether this could explain 
the worryingly small 0.5 mag offset between SOs and spirals. 

6 CONCLUSIONS 

We have determined the outer mass profiles in a sample of 18 

round galaxies, under the assumption that each is spherical 
and has constant anisotropy. The assumption of spherical 
symmetry is the more dangerous of the two: if a galaxy is 
fiattened along the line of sight, our spherical models under- 
estimate its mass, whereas varying tiic anisotropy iias little 
effect on the outer mass profile. Nevertheless, if we assume 
that a galaxy's anisotropy and the flattening of its light 
distribution and of its equipotential surfaces do not vary 
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Figure 3. Model fits for mass-to-light ratios, circular velocities and projected kinematics for the galaxies in our sample. The hatched 
regions give the 95% confidence bands for each quantity. For each galaxy an isotropic (/3 = 0) and a mildly radially anisotropic (,(3 = 0.3) 
model is constructed by varying the local mass-to-light ratio to fit the observed dispersion profile. The /i4 profiles plotted are model 
predictions, not fits. The PNe measurements used in our models of M32 and NGC 3379 are plotted as the crosses {v < 0) and open 
circles (i) > 0). 



strongly with radius then we may sensibly use our spherical 
models to check for evidence of a dark halo (i.e., a rise in 
mass-to- light ratio with radius). 

Our results are somewhat unsurprising. There is very 
clear evidence for massive dark haloes in largo ellipticals, 
whereas the the case for smaller galaxies is more ambiguous. 
Our early-type galaxies follow a TuUy-Fisher relation paral- 
lel to that for spiral galaxies, but offset by at least 0.8 mag 
(7-band). This is larger than the offset found by Neistein et 
al. (1999) for their sample of SO galaxies, but is reassuringly 
consistent with what one would expect from passive stellar 
evolution. 



The degeneracy betwoou aiiisotropy and shape pointed 
out in §2 (see also Magorriau 1999) means that we are un- 
able to place strong constraints on the anisotropies of our 
galaxies. It is often claimed that elliptical galaxies are radi- 
ally anisotropic, but there is surprisingly little solid evidence 
for this assertion. For example, detailed dynamical models 
have been constructed for the intermediate-luminosity ellip- 
ticals NGC 2434 (Rix et al. 1997), NGC 3379 (Gebhardt et 
al. 2000) and NGC 6703 (Gerhard et al. 1998). In all cases, 
the models were required to become radially anisotropic in 
order to fit the observed VP shapes, but in all cases it was 
unjustifiably assumed that the galaxy was either spherical 
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Figure 3. continued... 



or spheroidal - isotropic models with weak face-on disks can 
also provide reasonable fits to the VPs of these galaxies. 

Brighter ellipticals suffer less from this degeneracy, since 
they are unlikely to be disky. The strongest evidence for ra- 
dial anisotropy in an elliptical galaxy comes from Matthias 
& Gerhard's (1999) investigation of NGC 1600, a boxy E4 
galaxy. Their best-fit model has an anisotropy /3 ~ 0.5 in the 
equatorial plane and /? ~ along the symmetry axis. Since 
this galaxy is so strongly flattened in projection, it must bo 
close to edge on, minimizing the uncertainties in the shape of 
its intrinsic light distribution, and, therefore, its mass profile 
and anisotropy (provided one makes the reasonable assump- 
tion that its mass distribution is not significantly flatter than 
its light distribution). 

It is unlikely that one will ever be able to directly place 
strong constraints on the intrinsic shapes (and therefore 



anisotropics and masses) of round galaxies. To gain a proper 
understanding of the dynamical nature of elliptical galaxies, 
it would be more sensible to start with a thorough investi- 
gation of the shapes and profiles of the dark and luminous 
matter in fiattcned systems like NGC 1600. Models of these 
galaxies when "observed" face-on could then prove to be an 
interesting test of the nature of rounder galaxies. 
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Figure 3. continued. 
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APPENDIX A: THE DEGENERACY BETWEEN 
MASS, ANISOTROPY AND FLATTENING IN 
SCALE-FREE OBLATE GALAXY MODELS 

Gerhard (1993) has carried out an extensive investigation 
into how anisotropy affects the VPs of simple spherical 
galaxy models. Here we use scale-free models to compare 
the effects of anisotropy with the effects of flattening the 
galaxy's light distribution. For simplicity wc confine our at- 
tention to spheroidal models, and ignore the effects of disk- 
iness, boxiness or a radially varying ax:is ratio. 
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even stronger radial anisotropy can be mimicked by making 
the flattened galaxy slightly disky - the expressions in (A2) 
apply only for perfectly oblate galaxies. 

There are many galaxies whoso observed VPs arc more 
centrally peaked (larger hi) than the predictions from simple 
spherical isotropic models. One way of fitting these kinemat- 
ics is to construct spherical or spheroidal models with mild 
radial anisotropy (e.g., Rix et al. 1997; Gerhard et al. 1998; 
Gebhardt et al. 2000; Saglia et al. 2000). An alternative is 
to construct flattened models that are closer to isotropic, 
possibly with weak, face-on stellar disks. It is reasonable to 
invoke such disks to explain the kinematics of intermediate- 
luminosity ellipticals such as NGC 2434, NGC 3379 and 
NGC 6703, since many galaxies of similar luminosity are 
observed to be isotropic and disky when viewed edge on. A 
giant elliptical like NGC 1399, however, is much less likely 
to have a stellar disk, meaning that the degeneracy between 
anisotropy and shape is more strongly constrained to lie 
within the limits set by (A2) and (A3). 



Figure 5. Degeneracy between mass, anisotropy and flatten- 
ing in scale-free models. The panels on the left show how the 
Gauss-Hermite coefiicients of the VPs of face-on oblate isotropic 
models vary with axis ratio. The panels on the ligiit show how 
the VPs of spherical models depend on anisotropy. All models 
have a luminosity density profile j ~ r^** embedded in a spheri- 
cal isothermal halo, but similar results are found for other density 
profiles and for a Keplerian potential. The lines plot the crude fits 
described in the text. 

Our models have power-law density profiles 

2\ -"/2 
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with axis ratio q, and cither Kepler $(r) = —GM/r or loga- 
rithmic <3>(r) = Inr potentials. We construct two orthog- 
onal classes of models: spherical, anisotropic (g = 1, /? y^: 0) 
and flattened, isotropic (g 1, /3 = 0), viewed face-on. 

We calculate the VPs of the flattened, isotropic models 
using the method described in Magorrian & Binney (1994). 
For the spherical anisotropic models we solve equation (16) 
for fo{Q) and calculate the VPs by direct integration. This 
gives results that are in very good agreement with the 
moment-based method of §4. Figure 5 shows the Gauss- 
Hermite coefficients of the VPs of models with a = 4 in 
spherical isothermal halos with Vc = 1- Changing the axis 
ratio of one of the isotropic models by an amount Aq, the 
Gauss-Hermite cofflcicnts of its VPs change by 



^^0.8Ag 
cr 



A/l4 



-0.07Ag. 



(A2) 



On the other hand, on increasing a spherical galax:y's 
anisotropy by an amount A/3, its VPs change by 



^^-0.3A/3 

(T 

A/t4 ~ 0.12A/3. 



(A3) 



We have checked that these expressions hold for models with 
2 a < 4 in both Kepler and halo potentials. The effects of 



APPENDIX B: DISTINGUISHING BETWEEN 
MODELS WITH DARK HALOS AND MODELS 
WITH CONSTANT MASS-TO-LIGHT RATIOS 

Galaxies become very faint at large radii, meaning that 
absorption- line spectroscopy can be used to measure stel- 
lar velocity dispersions only out to some maximum radius 
ri --^ r-off. One way of probing a galaxy's mass distribution 
at radii larger than ri is to obtain radial velocities of a large 
number of PNe. Our aim in this appendix is to calculate 
the number of PNe that would be required to distinguish 
between constant-T and constant-i)c models at given con- 
fidence level. We assume that the anisotropy and intrinsic 
shape of each galaxy do not vary strongly with radius (§2). 

Let us write <7i for the galaxy's projected velocity dis- 
persion at ri and a2 for the projected luminosity-weighted 
velocity dispersion averaged over radii r > ri. The for- 
mer is estimated by absorption-line spectroscopy, the lat- 
ter from the PNe radial velocities. For n ~ Ves, we find 
that constant-T models of our galaxies predict a ratio TZ = 
d^jox — 0.7, whereas constant- models of course predict 
72- = 1. Wo write Mt for the model with constant mass-to- 
light ratio (7?. ~ 0.7) and Mhaio for the model with a halo 
(^=1). 

The available data D consist of an estimate a\ ± AtJi 
of (Ti, and radial velocity measurements v\ . . .vt^ of the N 
PNe with radii r > r\. We assume that \){a\\D) is Gaus- 
sian with mean a\ and dispersion Aai. The radial velocity 
measurements provide an estimate (t| = '^vf /N — (Aw)^ 
of (T2, where Aw is a typical error in the measured Vi. By 
the central limit theorem, when N is large p(a"2|-D) may be 
approximated by a Gaussian with mean 02 and dispersion 
AcT2 = o-2[(l + {Avla2f)l2Nf''^. Therefore 
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Setting a2 = itai , (72 = TZai and integrating over ai yields 

(B2) 

where fi = Affi I Ui are the fractional errors on the estimates 
of (Ti and (72. 

For an observed dispersion ratio Ti., Bayes' theorem 
gives the relative odds of Mhaio over Mr as 

p(Mhalo|7^) ^ p(7^|Mhalo) p(Mhalo) . 

p(MT|7e) p(7e|MT) p(Mt) ■ ^ ' 

Having no prior preference for either model, we take 
p(Mhaio) = p(Mt). Now suppose that Mr were the correct 
model. Then the fraction of observations that will correctly 
choose this model over Mhaio can be obtained by integrat- 
ing p(7?.|Mt) over all 72. for which the odds (B3) favour 
Mr- Performing this integral, we find that one must have 
(/i + llf^ ^ 0.1 if this fraction is to be larger than 95%. 
Therefore, distinguishing between constant-T and constant- 
tic models at the 95% confidence level requires first that 
one have an estimate of the velocity dispersion at ri ~ reft 
with fractional error /i < 10%. If this dispersion could be 
measured with 5% accuracy, then one would need to mea^ 
sure the mean dispersion outside ri with a fractional error 
/2 < 8.7%, which could be accomplished using iV ~ 70 PNe 
radial velocities. 
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